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Abstract. In this paper we develop a quantization method for flat compact manifolds 
based on path integrals. In this method the Hilbert space of holomorphic functions in the 
complexification of the manifold is used. This space is a reproducing kernel Hilbert space. 
A definition of the Feynman propagator, based on the reproducing property of this space, 
is proposed. In the K™ case the obtained results coincide with the known expressions. 
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1 Introduction 

The quantization of a system whose configuration space is a differentiable manifold is a far from 
exhausted problem. For instance, in the case of a Riemannian manifold, it is not known which 
quantization scheme best represents the curvature of the manifold, see [1, 2, 3, 17, 18, 19, 21, 
22, 30]. Such problems have both physical and mathematical motivations. On one side there is 
the problem of the existence of the structures involved in the quantization and, on the other, 
the possible applications to specific physical problems. 

When the manifold is a compact Lie group, Hall showed that the quantization of the group is 
naturally isomorphic to the quantization of the cotangent space of the group, i.e., a quantization 
on the phase space. The latter coincides with the complexification of the group, see [11, 14, 16]. 
In other words, there are two Hilbert space structures, one given by the functions on the group 
and the other by the functions in the complexification of the group, both structures being 
naturally related. These problems have been explored and generalized in several other works, 
see, e.g., [11, 12, 14, 16, 31]. 

It is known that the cotangent bundle of a Riemannian manifold admits a natural complex 
structure, compatible with the symplectic form and the natural lifting of the Riemannian metric 
if and only if the manifold is flat. In the Appendix A we show this result following the work of 
Gorbunov [6], see also [20]. 

In the case of an orientable connected compact flat Riemannian manifold (Euclidean space 
form) we show that there is a natural isomorphism between the Hilbert space of square integrable 
complex functions on the configuration space and the space of square integrable holomorphic 
functions on the phase space. The scalar products are defined with a measure given by the 
fundamental solution of the heat equation on each space. 

This space of holomorphic functions on the phase space turns out to be a reproducing kernel 
Hilbert space. Taking advantage of the existence of a reproducing kernel we obtain the above 
mentioned isomorphism and a path integral which coincides with the known expressions in the 
Euclidean case, see [5, 35]. 

In particular, the 3-dimensional orientable compact Euclidean space forms present a particu¬ 
lar interest for cosmology, since they could model the spatial part of the flat-universe models [4] . 
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See the most recent works of J. Levin et al., which seek to develop a plausible cosmological 
model using orientable compact Euclidean space forms of dimension 3 in agreement with results 
of observations made on the cosmic microwave background radiation [25, 26, 27, 28]. 

2 Flat Riemannian manifolds 

A connected complete flat Riemannian manifold is the quotient of M" by a subgroup F of the 
Euclidean group E{n), which has a free and properly discontinuous action. This is a part of 
a more general theorem by W. Killing and H. Hopf [34], 

Theorem 1. Let M be a Riemannian manifold of dimension n > 2 and zero curvature. Then M 
is complete and connected if and only if it is isometric to the quotient M n /T with T C E(n), 
where T acts freely and properly discontinuously. 

These manifolds are known as Euclidean space forms. The Euclidean group E(n) is the 
semidirect product of the groups 0(n ) and M n . An element 7 G T C E(n) is identified with 
7 = ( A , a), A G 0(n) and a G M n . The action of this group on an element x G M n is given by 
7 (x) = Ax + a. 

In one dimension the manifolds of this type are the real line M and the circle S 1 . In dimension 2 
there are five manifolds, the plane M 2 , the cylinder, the infinite Moebius strip, the torus and 
the Klein bottle. The torus and Klein bottle are both compact, while the torus is the only 
orientable one. In dimension 3 there are 18 types, 10 of which are compact, 6 orientable and 4 
non-orientable [23, 34]. In higher dimensions the number grows significantly; for example, there 
are 74 compact types in dimension 4. 

If T C E(n) is a lattice, then the quotient M n /T is called an n-torus, and is a compact 
Euclidean space form and a Lie group. 

It can be shown that every homogeneous Riemannian manifold is diffeomorphic to some Lie 
group but in general a space form is not necessarily homogeneous. In particular, when the space 
form is homogeneous we have the following theorem [34, p. 88 ]: 

Theorem 2. Let M be a connected homogeneous Riemannian manifold of dimension n and 
zero curvature, then it is isometric to the product M m x T n ~ m of a Euclidean space with a flat 
Riemannian torus. 

A discrete subgroup is a subgroup which is a discrete subset. If V is a closed subgroup of G, 
is called uniform if the quotient space G/T = {gT: g G G} is compact. 

The following theorem characterize discontinuous groups on Euclidean spaces. 

Theorem 3. Let T be a subgroup of the Euclidean group E{n). 

(i) T acts properly discontinuously on Euclidean space M n if, and only if, T is discrete on E{n). 
(ii) If T is closed, then it acts freely on M n if and only it is torsion free. 

{Hi) T acts properly discontinuously and with compact quotient on M n if, and only if, T is 
a discrete uniform subgroup of E(n). 

This paper focuses on orientable compact flat manifolds. 

Flat compact Riemannian manifolds of dimension n are quotients of polyhedra in by 
identifying faces (see [34, p. 99]). The interior of these polyhedra may be taken as a chart, 
which we call Q°. Functions defined on the manifold are functions on M n , which are invariant 
under the action of the group. 

An important invariant for a compact Euclidean space form is its volume. This can be defined 
in terms of a fundamental region for T in W 1 [29]. The volume of a space form M n /T is defined 
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to be the volume of any fundamental region. As a fundamental region for P we can take c 7 , the 
closure of the Dirichlet domain centered at 7 ( 0 ) 

c 7 := {x € M n ; || 7 ( 0 ) — x\\ < (^'(O) — x|| for every 7 ' € Tj, 

where 7 ( 0 ) is the action of 7 on 0 E M n . c 7 is an n-dimensional convex polyhedron in M n 
bounded by hyperplanes which are perpendicular bisectors of line segments [ 7 ( 0 ), 7 r ( 0 )]. Its 
boundary <9c 7 carries a locally finite decomposition into convex polyhedra of dimension n — 1. 
The space form M n /T is then obtained from by identifying points in <9c 7 which are equivalent 
modulo T. 

In particular, the six 3-dimensional orientable compact Euclidean space forms are the fol¬ 
lowing quotient spaces M 3 /Tj, i = (see [34, p. 117] and [23, p. 302]). The torus T 3 , 

which is constructed by identifying the opposite faces of a parallelepiped by translations, in this 
case Ti is generated by three translations t\ , t 2 , 1 3, in the direction of three linear independent 
vectors. Other four are obtained after gluing with a quarter turn, a half-turn, a one-sixth turn 
and a one-third turn. The last one is the Hantzsche-Wendt space which has a more complicated 
structure, 3 screw motions are needed. T 2 is generated by Ti and a screw motion a 2 = £ 3 , 
the faces of the translated parallelepiped are identified after a rotation of an angle of 77 T 3 is 
generated by Ti and a screw motion a 3 = £ 3 . r .4 is generated by Ti and a screw motion a 4 = £ 3 . 
T 5 is generated by Ti and a screw motion a 6 = t%. In the last one, the Hantzsche-Wendt space, 
T 6 results from T 2 by adding two further screw motions by an angle of tt. The manifolds M 3 /T 3 
and M 3 /T 5 are obtained from a lattice made by translating a hexagonal plane lattice a certain 
distance perpendicular to the plane and identifying opposite sides with the top rotated by 2ir/3 
and 7 t /3 respectively. 

The family {c 7 } forms a crystalline structure whose symmetry group contains T as a subgroup 
of finite index (see [34, p. 100]). We can choose the interior of Co (the cell corresponding to the 
identity of the group) as the chart Q° of M n /T. The crystalline structure can be generated by 
translation of a finite set of vectors defining the crystal lattice. This set forms a basis of M n . 
Dual basis vectors multiplied by 2 - 7 r are the basis of the reciprocal lattice, C. Let K be an 
element of the reciprocal lattice, then a function with the symmetry of this lattice has a Fourier 
expansion given by 

f(x) = c K e lK ' x - 
Kec 

This function is well defined on the manifold if it is also invariant under the action of T, i.e., 
(7 f)(x) = /( ix) = f{x) for all 7 G T. 


3 Hilbert spaces of quantization 

In a series of papers, Hall [11, 14, 16] showed that in a compact Lie group there is a natural 
isomorphism between the space of square integrable functions on the group, with a measure 
given by the fundamental solution of the heat equation at the identity, and the space of square 
integrable holomorphic functions in the complexification of the group (which is isomorphic to 
the cotangent space). Different but equivalent versions of this isomorphism can be found in the 
literature, see, e.g., [11, 12, 14, 16, 31]. 

For orientable compact Euclidean space forms we develop this isomorphism explicitly through 
integrals in Euclidean space and the solutions of the heat equations. Doing so paves the way for 
us to write an expression for the path integral in these spaces, which is the main result of this 
work. 
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3.1 Weighted heat kernel representation 


First we find the fundamental solution of the heat equation (heat kernel) in a compact flat Rie- 
mannian manifold. The solution of the heat equation defined on the manifold can be calculated 
by finding the solution on the chart given by the polyhedron Q°. Those solutions of the heat 
equation that are invariant under the action of the subgroup T are solutions on the manifold. 

If the manifold is M n , the heat equation 


dp?{x) 

dt 


\&Pt°(x), 


where A is the Laplacian, has the following fundamental solution 


( 1 ) 


Pt°(x) 


1 -(x-x 0 ) 2 /2t 

(27rf) n / 2 


( 2 ) 


This solution verifies the condition 


lim Pf°(x ) = 5(x — xq). (3) 

t-vo+ 

In the orientable compact case, this equation can be solved on the chart Q°, where we choose 
a point xq G Q°. We seek a function p^°(x) that verifies (1) with the condition (3). As was 
discussed above, the solution must be of the form 

p?W) = h E (4) 

K&C 


where we introduced the volume V of the cell because of normalization issues. Also, p^° (x) 
must be invariant under the action of T. Inserting the last expression in equation (1) and 
considering (3), we have 


P?{x) 


1 

V 


e iKix-x 0 )-K 2 t/2 

KeC 


(5) 


This expression is invariant under the group action. It is a consequence of the symmetry of the 
coefficients. Indeed, if 7 = (A, a) G T, then 

iK-'y(x—xo)—K 2 t/2 _ ’y ^ ^\K-A{x—xo)+U<-a-K 2 t/2 _ ^x 
^ KeC 


(7 pt) xo (x) = y e 
Kec 


The last equality results from the orthogonality of A and T having finite index in the symmetry 
group of the crystal. Then the fundamental solution in Q° is the fundamental solution in Q. 
We define a scalar product on Q by the following expression 

( f,9)Q= f f(x)g{x)pt°{x)dx. 

JQ 

This product gives a Hilbert space which we call L 2 (Q, p^°). We note that this product de¬ 
pends on the point xo, it is centered on a point of the manifold. The usual representation of 
quantum mechanics involves the space L 2 (Q) where the integration can be performed on the 
polyhedron Q° with the Lebesgue measure. The isometry between the two representations is 
given by 


/(*) 


fs{x) 


( 6 ) 


where fs is the corresponding function on L 2 (Q) and / on L 2 (Q, p^ 0 ). Here, we name the 
representation on L 2 (Q,p^°) weighted heat kernel representation. 
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3.2 Holomorphic representation 

The cotangent bundle of a flat manifold has a natural complex structure. See the Appendix A. 

We have seen that an appropriate chart for these manifolds is the interior of a polyhedron 
in W 1 , which we call Q°. Using this chart, the cotangent bundle has a natural chart Q° xR". The 
complex structure can be chosen so that the points of the polyhedron are the real coordinates 
of the complex manifold. 

Any real analytic function defined on Q° C has an analytic extension to Q ^ C C n , 
where Q ^ is Q = Q° xR". In particular, the fundamental solution (4) has an analytic extension 
on both variables x and xq. 

Consider the heat equation on the manifold 


dv;°(z) 

dt 




where z$ = xq + iyo is a point of Q° c and the Laplacian is taken with respect to the 2 n real 
variables x E Q° and y E of z = x + \y. Given the product structure of the chart, the 
fundamental solution can be calculated as the product of the solution on Q° by the solution 
on M n [7, 33]. 

In the Euclidean case, the heat kernel is given by 


«?(*) = 


(2vrt) r 


^-\z-z 0 \ 2 /2t. 


Furthermore, in the case of = Q° x M n described above we obtain 


= 


V (27Tt) rl / 2 


e -|Qf(^-Zo)| 2 /2t e iK-?R(z-z 0 )-K 2 t/2 

KeC 


(7) 


Then we define the scalar product of holomorphic functions 4>(z) and i/j(z) of as 

JQc 

The evaluation in t/2 has been done for convenience in order to obtain the isometry (11) below. 

We call T-LL 2 (Qc,^° 2 ) to the Hilbert space of square integrable holomorphic functions with 
this scalar product. In the Euclidean case (Q = M n ) it is the Segal-Bargmann space. 

The spaces of holomorphic functions shown in this paper are examples of reproducing ker¬ 
nel Hilbert spaces [13]. In these spaces there is a function I\(z,w), holomorphic in both ar¬ 
guments, i.e., K is holomorphic in z and antiholomorphic in w (holomorphic in w). For all 
4> £ T~LL 2 (Qr. t'f/g) the following identity is verified 


4>(z) = / K(z,w)4>(w)^? 2 (w)dw. (8) 

JQc 

In the Euclidean case the reproducing kernel can be obtained easily [13] 

K(z,w)= e* 0 /*. 

We consider linear operators on these spaces. Let A be represented by a kernel Ka(z,w ) as 
follows 



Ka{z, w)(t){w)iz^ 2 (w)(ivj. 


A<fi(z) 


(9) 
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Given an orthonormal basis of the Hilbert space {«i(£)}, * = 1, - -the reproducing kernel 
can be written 

OO 

K(z,w ) = y ^Uj(z)uj(w), (10) 

i =1 

and the corresponding kernel of the operator A is given by 

OO 

I\a(z,w) = y ^(Auj(z))uj(w). 

i =1 

Also, the composition of operators is associated with the following kernel 
I<ab(z, w ) = j K A (z, v)K b {v, w)v*° 2 (v) dv. 


3.3 Isometry between Hilbert spaces 

The spaces L 2 (Q,p x °) and T-LL 2 (Qc, v x ° 2 ) are naturally isomorphic. The isomorphism 

A t : L 2 (Q,tf°)^'HL 2 (Qc,v*° 2 ) (11) 

is given by integration over Q as follows 

Af{z)= f pt(x)f(x)dx, (12) 

JQ 


where pf(x) is the analytic extension of p*° (x) in the variable xq. 

The last expression is an isomorphism that can be explicitly tested as follows using the 
expressions for the kernels of the heat equations 


(A t f,A t g)Q c = f Atf{z)A t g{z)VtJ 2 (z)dz 
JQc 


/ f(x)g(x') / p^x)p z t (x')^°{z)dzdxdx' 

JQxQ JQc 

[ f{.x)g{x')p x t °(x')5(x - x')dxdx' = {f,g) Q . 

JQxQ 


The integral over Qc in the second line is evaluated using the expressions (5) and (7) along with 
the usual orthogonality relations. 

This form of the isometry is the analogous to what Hall calls Bt [11]. There are other forms, 
however we continue using (12) for the purpose of finding the path integral. 


3.4 The S 1 case 

Now we consider the case in which Q = S 1 . The fundamental solution p®°(6), centered on 9 o, 
satisfies the heat equation 


dp e t°{0) 

dt 


1 

2 


a S i p e m 


and converges to the Dirac delta 5(9 — 9q), for t —> 0 + . 
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The function p 00 is given by 

OO 

p e m= E Pk(t)e ik ^\ 

k= — OO 

and therefore the functions pk satisfy 
d Pk(t) 1, 2 m 

“dT = n{t) - 

and we finally obtain 






k( 8 — 6 o)—^k 2 t 


(13) 


fcez 


Then the scalar product of two functions / and g on S' 1 is 

OO 


u ’ s) = h 


E c m d n e i0(n - m) E e ik{d - 0o) ~ k t/2 d9 


m.n =—oo 


k=— oo 


= E c m d n e- i(m - n)0 °-( m - n)2t/2 = E E Cjdj-fc e - ik0 °- k2t/2 , 

m,n =—oo k=—o o \j=—oo ) 

where Cj and dj are the Fourier coefficients of / and respectively. 

The cotangent space is the cylinder S’ 1 xR. Thus, the complex manifold Sq can be represented 
by the chart (—vr, ir) x M viewed as the vertical strip in the complex plane. 

Holomorphic functions defined on the cylinder are holomorphic functions on C, which are 
also periodic on the coordinate corresponding to the real part. Using the Cauchy-Riemann 
conditions it is straightforward to see that they are of the form 

OO 

ifi{w) = E 

l = — OO 

Given that the heat kernel on a product manifold is the product of the respective heat 
kernels [7, 8], in this case in particular the heat kernel on S 1 x M is obtained from the heat 
kernel on S 1 (13) and the heat kernel on M (2), respectively. 

Then, the measure in this case is 


^°(z) = 


2-7T\/2irt 


a ( z ~ z o ) 2 

e 2t 


E 


3 in5R(z— zq)— n 2 t/2 


The Segal-Bargmann transform (12) of a function / can be obtained easily from his Fourier 
coefficients. It is given by 

1 P 7T 00 00 

tfi(z) = — / f(0) E e m{ - e -^~ n2t / 2 de = E c m e imz ~ m2t / 2 . 

~ n n =—00 m =—00 


Given that {e m:E /\/27r} with integer n is an orthonormal basis of L 2 (S’ 1 ), using (6) and (12) 

fin 

1 


we can obtain { <fi n }, which is an orthonormal basis of 77T 2 (S’ vn, 


Z = 


fTT QinXpZ^ 


dx. 


J —7T y/ Pt° ( X ) 

Then, by equation (10) the reproducing kernel is 

e ikx e ~iky 

fcez^ -7r V Pt°( x )Pt°(y ) 


k1z ' w)= ^LL. 


pt( x )pt(y) dx d v = ^ 


Pt( x )Pt(x) 
p? (z) 


d.x. 
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4 The Feynman path integral 

In this section we propose a definition of the Feynman integral in the holomorphic representation 
which is suitable for flat manifolds. Here, we take a heuristic motivation for our definition. 

First we consider the propagation of the wave function </> for a small time e. The propagation 
is governed by the Schrodinger equation with Hamiltonian operator H. We call Kh to the kernel 
of this Hamiltonian in the holomorphic representation. 

We obtain the evolution operator by exponentiating the Hamiltonian. Applying the evolution 
operator to (8), using (9) and introducing the series expansion of the exponential we obtain 



Thus, calling U e (j)(z) = e ieH (j>{z) and U e <j)(z) to the last integral, we have 
U e (j)(z) = U e cj>(z) + e^O, e). 

If we divide a time interval T into n equal parts, i.e., T = ne, we have 

U T = lim im., 


then 


U T <t>{z) = lim U^,J>{z) = lim U^, <f>(z) + - F(z,T/n) 

n—>• oo / n—>• oo \ ' \ Tl J 


TV 


= U T (f>(z). 


We define a Feynman propagator for a finite time T by dividing the interval [0, T] into n equal 
subintervals and taking the limit n —> oo. Then, it follows 


-i£ E h { Zj - Zj .o AT 

G(zt, z 0 ; T) = Jhm^ / ^ e K(zj,Zj- 1 ) v^ 2 (zj)dzj, 

,wo ° IQI 1 j=1 j=1 


(14) 


where we introduce the normal symbol of the Hamiltonian 

,, — \ K h (z,w) 

h{z, w ) = ~rr ( — it . 

K (z, w) 

The above expression can also be calculated in the Euclidean case. Here, we take the Hamiltonian 

h=4, 

corresponding to a renormalized one-dimensional harmonic oscillator. The kernel of this Hamil¬ 
tonian is Kh(z,w) = zwe zw and the normal symbol h(z,w) = zw. 

The Feynman propagator (14) in this case is 


G(z t ,zo]T) = lim / &- 1 

n—>OC J^n-i 


n n — 1 1 

E zjZj- i-i£ E Zjtj- 1- E M 2 d Zj 


j=l J=1 yxz 3 


3= 1 


2t r 
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Then, by regrouping terms, we have 


G(z t ,z 0 -,T) = lim 


n —>og 


J 



where again e = T/n. 

This expression is the known Feynman integral 



shown in [5, 35]. 

5 Discussion 

The cotangent bundle of a Riemannian manifold admits a complex structure compatible with 


the symplectic form an the natural lifting of the Riemannian metric if and only if the manifold 


is flat, i.e., this structure is not integrable unless the curvature is zero, see Appendix A. 

The restriction on the integrability of the almost-complex structure can be circumvented 
by changing the metric. In [6], in the context of deformation quantization, Gorbunov et al. 
constructed a formally Kahler metric on the cotangent space by adding powers of momentum, 
obtaining an integrable Kahler structure. 

Furthermore, it has been shown that a complex manifold can be constructed in a neighbor¬ 
hood of the null section of the tangent bundle of a real Riemannian manifold. These are called 
Grauert tubes [9, 10, 24, 32], The complex structures defined therein, called adapted complex 
structures, are compatible with the symplectic structure leading to a Kahler manifold [15]. In 
certain cases this complex structure exists throughout the tangent, for example when the base 
manifold is a compact Lie group with a bi-invariant metric. 

Our paper is focused on flat Riemannian manifolds, specifically on Euclidean space forms. 
They present a particular interest in cosmology because they could model the spatial part of 
flat universe models, see [4, 25]. The formulas presented in this paper are applicable only to flat 
manifolds. However the method that leads to the formulation of the path integral is interesting 
in itself. Moreover, it allows exploring the quantization of space-forms, a topic sparsely discussed 
in the literature. Finally, the method presented here could be useful for generalizations in future 
works. 

A Complex structure in phase space 

It is known that the cotangent bundle of a Riemannian manifold admits a natural complex 
structure, compatible with the symplectic form and the natural lifting of the Riemannian metric 
if and only if the manifold is flat [6, 20]. 

First we consider the general case. Let Q be a finite-dimensional orientable real-analytic 
Riemannian manifold and let T*Q be the cotangent bundle with projection it: T*Q —> Q and 
canonical symplectic form ui. 

The metric in Q can be naturally lifted to the cotangent space T*Q as follows. Let a(t) = 
(qi(t),pi(t)) and /3(f) = (q 2 (t),P 2 (t)) be curves on T*Q such that 


a(0) = /3(0) = (g,p) = m, 


a'( 0) = V, and /3'(0) = W, 
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i.e., V and W are tangent vectors on m. We denote by a the metric on Q, and by a$ the 
isomorphism induced by a between T*Q and TQ. Then we have a metric G on T*Q given by 
the following expression 

Gm(V,W) = a q (TttV, TttW) + (^^(0), a“^(0)) . 

Here, T7r is the tangent application of the projection and D represents the covariant derivative. 

Thus, we have two different structures on T*Q\ the canonical symplectic structure ui and the 
Riemannian metric G. A third one appears naturally. It is the almost-complex structure J. If 
ufi \ T*T*Q — > TT*Q is the isomorphism induced by ui, and G^: TT*Q -» T*T*Q the isomor¬ 
phism induced by G , then J is given by J = uj^G". It means that given the fields V and W, 
they verify 


G(V,W) =u(V,JW). 

At each point m the map J m : T m T*Q —> T m T*Q verifies J^V = — V. J is compatible with u 
and G. (J, G, u) form what is known as a compatible triple. 

We are now interested in finding appropriate local complex coordinates. The complexified 
tangent bundle T C T*Q splits as follows 

T C T*Q = T (1 ’ 0) T*Q © T ( -°’ 1 ' ) T*Q. 


Here T^ 1,0 ^T*Q and T^ 0,1 ^T*Q are the images of the projections n + and n given by 


n ± 


1 © iJ 
2 


If we choose a vector V = (g 1 ,..., q n ,pi ,... ,p n ) E T m T*Q , then n + is a natural isomorphism 
between the real tangent space T m T*Q and the holomorphic tangent space Tm’°^T*Q. Explicitly, 
we have 


U + V 



where z l = q l + icr* m (p m - p k T k ml q l ). 

The corresponding holomorphic vector fields are 


d _ 1 ( d - 9 \ 

dz l 2 \ dq l ^ k dp> ^ dpi J 


where ©j and T ? k are the matrix coefficients of the metric and the Christoffel symbols respec¬ 
tively. Henceforth, the Einstein summation convention is used. 

If we take the Lie bracket of the above-mentioned fields we obtain 


Then, by the Niremberg-Newlander theorem [20], the distribution is integrable if and only if 
the curvature tensor of the metric a is identically null, see Gorbunov et al. [6] . 



d_ d_ 
~dzi' dzi 


= iR% jPm a lk I 


d _d_ 
d¥ 
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